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Controlling the non-Markovian dynamics of open quantum systems is essential in quantum in¬ 
formation technology since it plays a crucial role in preserving quantum memory. Albeit in many 
realistic scenarios the quantum system can simultaneously interact with composite environments, 
this condition remains little understood, particularly regarding the effect of the coupling between en¬ 
vironmental parts. We analyze the non-Markovian behavior of a qubit interacting at the same time 
with two coupled single-mode cavities which in turn dissipate into memoryless or memory-keeping 
reservoirs. We show that increasing the control parameter, that is the two-mode coupling, allows for 
triggering and enhancing a non-Markovian dynamics for the qubit starting from a Markovian one in 
absence of coupling. Surprisingly, if the qubit dynamics is non-Markovian for zero control param¬ 
eter, increasing the latter enables multiple transitions from non-Markovian to Markovian regimes. 
These results hold independently on the nature of the reservoirs. This work highlights that suit¬ 
ably engineering the coupling between parts of a compound environment can efficiently harness the 
quantum memory, stored in a qubit, based on non-Markovianity. 

PACS numbers: 03.65.Yz, 03.67.-a 


I. INTRODUCTION 

A thorough understanding of the dynamics of an open 
quantum system has experienced a long term pursuit 
[1] and nowadays attracts ever-increasing attention due 
to the development of quantum information technology 
HIS] that employs open quantum systems as basic re¬ 
source. In the theory of open quantum systems, the 
non-Markovian dynamics is one of main concerns being 
linked to the preservation of quantum memory stored in 
a quantum system It arises in many realistic 

situations HES], and also proves useful in quantum in¬ 
formation processing such as quantum state engineering, 
quantum channel capacity and quantum control 
[28] . The degree of a non-Markovian evolution, the so- 
called non-Markovianity, can be quantified by different 
measures based on dynamical features of the system ca¬ 
pable to grasp the memory effects of the environment 
on the system evolution [28ll35] . So far, many factors 
that can trigger and modify the non-Markovian dynam¬ 
ics have been found, as strong system-environment cou¬ 
pling, structured reservoirs, low temperatures, and initial 
system-environment correlations [Dimisnini- Apart 
from these mechanisms, some other peculiar conditions 
such as classical environments miiig and environmental 
initial correlations [43l [44] have also been predicted and 
experimentally demonstrated [laiii] enabling emergence 
of non-Markovianity. 

In the conventional study, one usually considers the 
quantum system being coupled to a single environment 
mills]. However, in several realistic scenarios the system 


may be simultaneously influenced by many environments 
[4SH49]. For instance, in a quantum dot the electron spin 
may be affected strongly by the surrounding nuclei and 
weakly by the phonons [45] . The neighbor nitrogen impu¬ 
rities constitute the principal bath for a nitrogen-vacancy 
center, while the carbon-13 nuclear spins also have some 
interaction with it [46] . A similar situation also occurs for 
a single-donor electron spin in silicon [47]. Motivated by 
these practical situations [45H47] , some efforts have been 
devoted to study the effects of multiple environments on 
the dynamics of an open system [48H52] . Quantum inter¬ 
ference effects have been found to occur between indepen¬ 
dent reservoirs when all of them interact with a quantum 
system and are in non-Markovian regimes, which quali¬ 
tatively modify the dynamics of the interested system 
[48] . The dynamics of a spin simultaneously coupled to 
two decoherence channels, one Markovian and the other 
non-Markovian, has been analyzed with respect to the 
different decoherence mechanisms [49]. As is known, a 
qubit (i.e., a two-level system) interacting with a single 
vacuum bosonic reservoir may exhibit Markovian or non- 
Markovian dynamics depending on the strength of the 
system-reservoir coupling [T]. By contrast, if the qubit 
simultaneously interacts with several reservoirs, its dy¬ 
namics can be always non-Markovian provided that the 
number of the contributing reservoirs is greater than a 
critical value [50]. The dynamics of a qubit coupled to 
a hierarchical environment made of a single-mode cavity 
and a structured reservoir with Lorentzian spectral den¬ 
sity has been studied, showing that a shorter (longer) 
memory time of the reservoir does not generally mean 
a smaller (larger) non-Markovianity of the system [51]. 
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A new analytical method based on a phase space rep¬ 
resentation of the system density matrix has been also 
proposed to study the dynamics of a discrete system in 
a two-tiered non-Markovian environment [52] . 

In the treatments of composite environments [481150] . 
the role of the coupling between environmental parts is 
not typically taken into account. Despite this, the knowl¬ 
edge of how this environmental parameter influences the 
non-Markovian character of an open quantum system 
would provide insightful developments for engineering 
and controlling quantum memories for possible techno¬ 
logical applications. Therefore, this aspect deserves care¬ 
ful investigation, possibly starting from a paradigmatic 
model where it can simply emerge and be understood. 
Here we choose a model which complies with this re¬ 
quirement, namely a qubit interacting at the same time 
with two coupled single-mode cavities which in turn dissi¬ 
pate photons into their own memory less (Markovian) or 
memory-keeping (non-Markovian) reservoirs. This sys¬ 
tem finds its natural implementation in nowadays tech¬ 
nologies of circuit quantum electrodynamics [53| and also 
in simulating all-optical setups [54|. We shall show that 
the coupling strength between the two modes can harness 
the qubit non-Markovianity in different and even counter¬ 
intuitive ways, independently of the nature of the reser¬ 
voirs. The paper is organized as follows. In Sec. |TI| we 
describe the system of interest. In Sec. Ill we study the 
case of memory less reservoirs, while in Sec. |IV| we inves¬ 
tigate the case of memory-keeping reservoirs. In Sec. [V| 
we summarize our conclusions. 


II. THE SYSTEM 

Our global system consists of a qubit s simultaneously 
interacting with two environments £i and £ 2 - To fix the 
ideas, we model each environment £n {n = 1,2) as a 
bosonic mode rrin decaying to a zero-temperature bosonic 
reservoir 7^^, as depicted in Fig. The qubit is here 
meant as the quantum memory whose efficiency is to be 
quantified by its non-Markovianity, that is by the degree 
of non-Markovian evolution. The interaction of the two 
environments is due to the coupling of the two bosonic 
cavity modes, which instead plays the role of a control 
parameter for the non-Markovianity of the qubit. The 
coupling strength of the qubit with each mode rrin is 
while Q denotes the coupling between the two modes. 
For the sake of simplicity, we assume that the two modes 
have the same frequency uJc which in turn is equal to the 
qubit transition frequency ujq, that is ujq = cJc- 
The total Hamiltonian is given hy {h = 1) 
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n=l 

"b-Hmim2 5 (1) 

where Hs = {uJo/2)az is the qubit Hamiltonian, Hrun = 
LOcahan the mode Hamiltonian, Hn„ = ^n,kb\ )})n,k 



FIG. 1: Pictorial description of the global system. A 
qubit, i.e. a two-level system (TLS), simultaneously interacts 
with two environments £n (n = 1,2), each containing a single 
cavity mode rrin that in turns is affected by a vacuum (zero 
temperature) reservoir IZn- The qubit is directly coupled to 
each mode rrin with strength t^n- The two cavity modes mi, 
m 2 are coupled with strength £t. 


the reservoir Hamiltonian, Hsmn = Rn{^h^- + 
the qubit-mode interaction Hamiltonian, — 

^n,/c(^n^n/c + ^n^n,/c) ftie modc-rescrvoir interaction 

Hamiltonian and Hmim 2 = f7(d|d2 + did^) the interac¬ 
tion Hamiltonian between the two modes. In the expres¬ 
sions above az = |1) (1| — |0) (0| is a Pauli operator for 
the system with transition frequency ujq, (J± represent 
the raising and lowering operators of the qubit, (d^) 
the annihilation (creation) operator of mode rUn- Fur¬ 
thermore, in the Hamiltonians involving the reservoirs 
bn,k (bl^ k) annihilation (creation) operator of field 

mode k with frequency ujri,k of reservoir 7^^, and g^^k de¬ 
notes the coupling of the mode rrin with the mode k of 
its own reservoir 7^^. In the interaction picture, the total 
Hamiltonian can be expressed as 

2 

-^int ~ ^ ^ Hgnriri T Hni-^m2 

n=l 

2 

n=l k 

where An,k = ^n,k - ^o- 

The reservoirs IZn of the global system can be ei¬ 
ther memoryless (Markovian) or memory-keeping (non- 
Markovian). Depending on the kind of reservoir, different 
methods are used to obtain the reduced dynamics of the 
qubit. In the following two sections we study these two 
cases. 
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III. MEMORYLESS RESERVOIRS 


In this section we consider both the reservoirs IZn as 
vacuum Markovian ones, their correlation times being 
much smaller than the single-mode relaxation times . Al - 
though our system can be exactly solved (see Sec. IV), 
we first treat it under the Markov approximation since 
this analysis constitutes a strategical first step in order to 
strongly evidence the crucial role of the two-mode cou¬ 
pling parameter to harness quantum non-Markovianity 
for the dynamics of the qubit even under this condition. 
In this case, the density operator p(t) of the qubit plus 
the two modes obeys the following master equation [T] 


p{t) = -i[H,p{t)] 

^ ^ r 

- X] Y “ 2anp{t)al + p{t)alan], (3) 

n=l 

where p{t) = dp{t)/dt, H is given by Eq. Q with¬ 
out the terms involving the reservoirs and denotes 
the decay rate of the mode rrin- We initially take the 
qubit in its excited state |1)^ and both modes in the 
ground states |00)^^ so that the initial overall state 
is p(0) = 1100) (100|, where the first, second and third 
element correspond to the qubit 5, mode mi and mode 
m 2 , respectively. Since there exist at most one excitation 
in the total system at any time, we can make the ansatz 
for p{t) in the form 


p{t) = (1 - Ht)) IV’(i)) (V’(i)l + A(i) |000) (000|, (4) 


where 0 < A(t) < 1 with A(0) = 0 and |'0(t)) = 
h{t) 1100) -t- ci{t) |010) -|- C 2 {t) |001) with h{0) = 1 and 
ci(0) = C2(0) = 0. It is convenient to introduce the un¬ 
normalized state vector m 


(a) 




n/r 


FIG. 2: Non-Markovianity measure A/^ as a function of the 
coupling constant Q/F between the two modes for (a) weak 
and (b) strong system-mode coupling regimes. 


^W) = \/l - A(t) \ip{t)) 

= ft(t)|100)+Si(i)|010)+2-2(^)1001), (5) 

where h{t) = — X{t)h{t) represents the probability 

amplitude of the qubit and c^(t) = — X{t)cn{t) that 

of the mode mn being in their excited states. In terms 
of the unnormalized state vector we then get 

p(t) = |^(t)){^(t)|+A(t)|000)(000|. (6) 

Inserting this expression in Eq. the time-dependent 
amplitudes h(t), ci(t), 02 ( 1 ) of Eq. (O are determined by 
a set of differential equations as 

ujQh{t) -h nici{t) + n2C2{t), 

Ci{t) + K^{t) + 0C2(1), 

- ^r2^ C2(i) + K2^(1) + eci(l). (7) 


The above differential equations can be solved by means 
of standard Laplace transformations combined with nu¬ 
merical simulations to obtain the reduced density opera¬ 
tors of the qubit as well as of each of the modes. 

To quantify the non-Markovianity we adopt a measure 
based on the dynamics of the trace distance between two 
different initial states pi(0) and ^ 2 ( 0 ) of an open system 
[ 29 ]. A Markovian evolution can never increase the dis- 
tinguishability between different initial states in terms of 
their trace distance, hence a nonmonotonic time behav¬ 
ior of the latter would signify non-Markovian dynamics of 
the system. Such a measure is consistent with the inter¬ 
pretation of non-Markovianity in terms of a backflow of 
information from the environment to the system, which 
is responsible for the distance (state-distinguishability) 
growth. Based on this concept, the non-Markovianity 
can be quantified by a measure M defined as [29] 

J\f = max / cr[t,pi( 0 ),p 2 ( 0 )]dt, (8) 

Pl(0),P2(0) Jcr>0 

in which cr[t, pi(0), p2(0)] = dD[pi(t)^ p 2 {t)]/dt is the rate 


.dhit) 

i —;— 

dt 

. dci it) 

% -;- 

dt 

. dC2 {t) 
dt 
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of change of the trace distance given by 

D[pi{t),p2{t)] = ^Ti\pi{t) - P2{t)\, (9) 

where |^| = VA''iA. In order to evaluate J\f, we have to 
find a specific pair of optimal initial states to maximize 
the time derivative of the trace distance. In Ref. [56] . 
it is proved that the pair of optimal states is associated 
with two antipodal pure states on the surface of the Bloch 
sphere. We thus adopt a pair of initial states pi, 2 ( 0 ) = 
|■*/’l. 2 ( 0 )) (■i/’i. 2 ( 0)1 with |V’i, 2 ( 0 )) = (| 0 ) ± |l))\b as the 
optimal ones throughout the paper. This allows us to 
obtain the time derivative of th^ trace distance in the 
simple form cr[t, pi( 0 ), ^ 2 ( 6 )] = d\h{t)\/dt. 

In the absence of the coupling between mi and m 2 , 
the qubit exhibits Markovian dynamics when the cou¬ 
plings of the qubit with the two modes in terms of 
Ki/Ti and K 2/'^2 are weak. In this case, we show that 
the introduction of mode-mode coupling with sufficient 
strength Q can transform the Markovian dynamics to 
the non-Markovian one. In Fig. |^a), we plot the non- 
Markovianity A/" as a function of the scaled coupling 
strength f^/F between the two modes for different val¬ 
ues of hii and K 2 (ri=r 2 =r is assumed along the paper). 
As shown in the figure, the system exhibits Markovian 
dynamics, individuated by A/" = 0 , until the two modes 
are weakly coupled and below a certain threshold. How¬ 
ever, when Q/r exceeds this threshold the Markovian dy¬ 
namics of the qubit changes to non-Markovian one (i.e., 
Af > 0). In general, non-Markovianity increases with 
r^/F for the given values of ki and 1^2 and is also propor¬ 
tional to /^i, tv 2 foi* a fixed f^/F. Therefore, the coupling 
of the two modes can trigger the non-Markovian dynam¬ 
ics of the system. 

On the other hand, if the qubit-mode couplings /^i, k ,2 
are strong the qubit exhibits non-Markovian dynamics 
without the need of mode-mode coupling. Under these 
conditions, how the additional coupling of the two modes 
influences the system non-Markovianity is to be revealed. 
From the above discussion about the case of weak qubit¬ 
mode couplings, one might expect that the mode-mode 
coupling would enhance the non-Markovianity of the sys¬ 
tem. However, as shown in Fig. |^b), for different /^i, k ,2 
the relation between the non-Markovianity Af and the 
mode-mode coupling strength is nonmonotonic. Increas¬ 
ing U/F from zero, the non-Markovianity first diminishes 
to a minimal value and then rises. Remarkably, for some 
smaller values of ki and /^ 2 , the non-Markovianity can 
even decrease to zero (Markovian regime), remain zero 
for a finite range of O/F and then recover nonzero values 
with a further increase of U/F. The mode-mode coupling 
is thus able not only to enhance the memory effects of 
the overall environment but also to restrain them. 

For a comprehensive understanding of the effects of 
system-mode and mode-mode couplings on the qubit dy¬ 
namics, in Fig. we display the phase diagram in the hz-fl 
plane {ki = K 2 = k,) of the transitions between Marko¬ 
vian and non-Markovian dynamics. In the strong qubit¬ 



FIG. 3: Phase diagram in the hi/T-Q/T plane for the crossover 
between Markovian and non-Markovian dynamics. The col¬ 
ored regions represent the non-Markovian dynamics, while the 
white regions denote the Markovian dynamics. The dotted 
line near k/F = 0.18 divides the weak and strong couplings 
between the system and the modes and above (below) which 
is the strong (weak) system-modes coupling regime. 


mode coupling regime (above the dotted line), the qubit 
may experience two transitions: from non-Markovian to 
Markovian and again to non-Markovian dynamics (e.g., 
from point A to B and then to C) with an increase of 
O. In the weak qubit-mode coupling regime (below the 
dotted line), an increase of the mode-mode coupling can 
drive the Markovian dynamics to the non-Markovian one 
(e.g., from point D to E). Moreover, the smaller the /d, 
the larger the O required to activate the non-Markovian 
dynamics. 

So far, we have shown that the mode-mode coupling 
can trigger and modify the non-Markovianity of the sys¬ 
tem. In fact, the two modes can be regarded as effective 
memories of the overall environment since the compen¬ 
sated rate W{t) of their population changes, given by 

m 


W{t) 


d{EL, |Cn(i)p) 

dt 


n=l 


( 10 ) 


completely determines the qubit non-Markovian dynam¬ 
ics. The meaning of Eq. (10) can be explained as follows. 
The energy dissipations are one-way from the modes to 
their memoryless reservoirs so that Fi, r 2 are always 
positive. If the energy of the two modes decreases (i.e.. 


d /dt < 0 ) and this decrease is not com¬ 

pensated by the dissipation of the modes, quantified by 
Y/n=i rn|cn(^)P, then W{t) < 0. This situation can only 
happen when part of the two-mode energy has come back 
to the quantum system, thus a negative value of W{t) 
identifies a back-action (or back-flow of information) and 
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FIG. 4: Dynamics of the trace distance D[pi(t)^ p2{t)] (top plane) and the witness W(t) (unit s of Eq.(lO) (bottom plane), 
for — O.Sr and D = 0 (a), D = F (b) and Q = 2F (c). 


therefore onset of non-Markovianity. Therefore, the re¬ 
sults discussed above indicate that the coupling between 
the two environmental modes, the memories, can enhance 
or inhibit their own memory effects on the qubit. As a 
further verification, in Fig. we compare the dynamics 
of the trace distance D[pi{t)^ p2{t)] of Eq. ^ and of the 
rate W (t) to directly demonstrate that the negativity of 
W{t) assesses qubit non-Markovianity. To this aim, we 
choose three points in the k-Q phase diagram (Fig. 
with the same k = 0.3F, while = 0, F and 2F, respec¬ 
tively. These points pass through the regime transitions 
non-Markovian ^ Markovian ^ non-Markovian as O in¬ 
creases. As expected, D[pi{t), p2(j^] exhibit oscillations 
for = 0 and = 2F (see Fig. panels (a) and (c)), 
while asymptotically decay to zero for = F (Fig.|^b)). 
The witness W{t) becomes negative at the points where 
D[pi{t), p2{t)] begins to grow and remains negative dur¬ 
ing the whole time interval when D[pi{t), p2{t)] increases, 
which entails an information backflow from the modes 
to the qubit. Differently, W{t) remains positive when 
D[pi{t), P2{t)] asymptotically decays, as in Fig. [^b). It 
is worth to notice that, although the qubit undergoes 
non-Markovian dynamics for both 0 = 0 (Fig. |^a)) and 
O = 2F (Fig. Kc)), the dynamical curves of the trace 
distance D[pi{tJ^ P2{t)] (state distinguishability) are very 
different regarding the points when it starts increasing, 
implying different mechanisms of information backflows 
in the two cases. 


IV. MEMORY-KEEPING RESERVOIRS 

In the above section, we have considered a qubit s in¬ 
teracting with two coupling modes mi, m 2 which are 
dissipated respectively by two memoryless reservoirs IZi , 
7^2- Under these conditions, we have seen that the two 


modes are fully responsible for the memory effects of the 
overall environment on the qubit and their coupling can 
modify this effect. However, if the two modes are only 
components of the compound memory of the overall envi¬ 
ronment, the way their coupling changes the overall mem¬ 
ory effects on the qubit is to be explored. We accomplish 
this analysis in this section. To this purpose, we consider 
a more complex situation where the coupling modes are 
dissipated by structured reservoirs T^i, 7^2 exhibiting in¬ 
herent memory effects and are therefore non-Markovian 

BE]. 

We again assume the two modes and the correspond¬ 
ing reservoirs are initially in their ground states with 
only up to one excitation present in the total system. 
The qubit can initially be in a (normalized) superposi¬ 
tion \'ip)s = co(0)|0)s + /i(0)|l)s, so that the initial total 
state is |$(0)) = |tA)s|00)mim2 |00)7^J7^2 with |o).^^ = 
Ylk The total evolved pure state then reads 

\m) = [Co(0)|0), + Mi)|l)j|00)„^„j00)^^,,^ 

+ |0)s [Cl W |00 )k^7J2 

2 

+ ( 11 ) 

n=l k 

where |l/c)7^^ = |0 • • • 1/. • • • 0)7^^ means that there is one 
excitation in the kth mode of the reservoir IZn and n is 
the complementary of n (i.e., n = 2 if n = 1 and vicev- 
ersa). The initial conditions of the coefficients appearing 
in |i>(t)) are ci(0) = C2(0) = Cn,/c(0) = 0. From the 
Schrodinger equation [T], the time evolution of the total 
system in the interaction picture with the Hamiltonian of 
Eq. <§ is determined by the following differential equa¬ 
tions 

h{t) = - iK2C2{t), 
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FIG. 5 : The non-Markovianity A/^ as a function of the cou¬ 
pling constant Q/7 between two modes that are dissipated by 
Lorentzian reservoirs with 71 =72 =7. (a) and (b) demon¬ 
strate weak and strong couplings between the qubit and the 
modes in terms of and 1^2 ■ The other parameters are chosen 
as Ai = A2 = 0.57 in (a) and Ai = A2 = O.87 in (b). 


Ci(i) = -iK,ih{t)-iQ,C2{t) 

C2{t) = -iK2h{t) - iflciit) - 

C2,k{t) = (12) 

Integrating the last two equations with the initial condi¬ 
tion Cn,/c(0) = 0 (n = 1,2) and inserting their solutions 
into the second and third equation above, one obtains two 
integro-differential equations for the amplitudes ci (t) and 
C2{t) 


JO u 


The sum above equa¬ 

tions is recognized as the correlation function /^(t — t') 
of the reservoir 7^^, which in the limit of a large number 
of modes can be changed into an integration in terms of 
the spectral density Jn{^) as [T] 

fn{t-t') = j dioJn{u})exp[i{LOo-io){t-t')]. ( 14 ) 


We take each reservoir IZn with a Lorentzian spectral 
density Jn{uj) = 7nAy{27r[(u; - cjq)^ + A^]}, where 7^ 
is the mode-reservoir coupling strength and A“^ the 
reservoir correlation time [H [5]. The two-point corre¬ 
lation function of Eq. ( [l4| ) can be then expressed as 
/n(^) = ^7nAn exp(— Anf^. Therefore, the amplitudes 
h{t), ci(t) and C2(t) can be obtained by solving the first 
one of Eqs. together with Eq. ( p^ by using the 

standard Laplace transform technique. The reduced dy¬ 
namics of the qubit and the of other parts of the overall 
system are then determined by tracing out the opportune 
degrees of freedom from the evolved total state |^(t)) of 


Eq. 0. 

In Eig. we plot the non-Mar kovianity A/* as a func¬ 
tion of the scaled mode-mode coupling strength fl/y for 
different values of system-mode couplings /^i, k ,2 and 
assuming 71=72=7* As shown in Eig. |^a), when the 
memory effects of the two reservoirs alone (that is with 
fl/j = 0) are not sufficient to make the system ex¬ 
perience non-Markovian dynamics, the introduction of 
mode-mode coupling with sufficient strength can drive 
the Markovian dynamics to the non-Markovian one. The 
non-Mar kovianity is moreover proportional to the cou¬ 
pling strength for given values of k,^ and A^. This im¬ 
plies that the coupling of the modes, as constituents of 
the compound memory of the overall environment, can 
further enhance the memory effects of the latter on the 
qubit. On the other hand, if the system already under¬ 
goes non-Markovian dynamics without mode-mode cou¬ 
pling, the relation between the non-Mar kovianity and the 
mode-mode coupling results to be nonmonotonic. The 
mode-mode coupling reduces the non-Markovianity and 
can even transform a non-Markovian dynamics {Af > 0) 
into a Markovian one {Af = 0). Nevertheless, the fur¬ 
ther increase of the coupling strength Q can recover and 
increase the non-Markovianity, as shown in Eig. ib). 
These behaviors are analogous to the ones found before 
for the case of memory less reservoirs, a slight difference 
being that here the non-Markovianity exhibits oscilla¬ 
tions as Cl/j increases before reaching its minimal value 
(see Eigs. [^b) and|^b)). 


Ci{t) 


C2(t) 


—iKih{t) — iClc2{t) 


- [ 

•^0 k 

—iK2h{t) — iOci(t) 


2 -iAi,kit-t') 




V. CONCLUSION 

In this paper we have addressed the study of the effects 
of the coupling between two parts of a multiple environ- 
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ment on the dynamics of a quantum system. In par¬ 
ticular, we have considered a qubit (the quantum mem¬ 
ory) simultaneously interacting with two coupled bosonic 
modes (the control devices) which are in turns dissipated 
into memory less or memory-keeping reservoirs. In the 
case of memory less reservoirs, we have proven that the 
two cavity modes play the role of unique memory sources 
of the overall environment for the qubit and their cou¬ 
pling can be thus viewed as a coupling between two quan¬ 
tum memory sources. In the case of memory-keeping 
reservoirs, the two modes are instead the constituent 
parts of the total memory source of the overall environ¬ 
ment and their coupling can be now meant as a coupling 
between two partial memory sources for the qubit. We 
have shown that in both cases the Markovian dynamics 
of the qubit, existing without the mode-mode coupling 
(the control parameter), can become non-Markovian by 
adjusting the control parameter over a certain threshold. 
Moreover, higher values of the control parameter enable 
larger non-Markovianity for the qubit. Differently, when 
the qubit evolution is already in a non-Markovian regime 
for a zero mode-mode coupling, a nonmonotonic relation¬ 
ship arises between non-Markovianity and control param¬ 
eter. Namely, multiple crossovers from non-Markovian to 
Markovian regimes may occur by increasing the mode¬ 
mode coupling. This may appear surprising since, on 
the basis that increasing the coupling between memory 
sources for the qubit entails a transition from Markovian 
to non-Markovian regimes for the qubit, one expects that 
an increasing of the control parameter always induces an 
enhancement of memory effects on the qubit dynamics. 

We remark that the behaviors above happen indepen¬ 
dently of the nature of the reservoirs. Our findings ev¬ 
idence that when the environment is composite the un¬ 


derlying physical mechanisms may be counterintuitive. 
The environmental coupling thus reveals as a powerful 
and effective tool to activate and harness quantum non- 
Markovianity of open systems. It is worth to notice that 
our system has the advantage to make it emerge in a 
clear way the effects of this coupling on the dynamics 
of a quantum system and, at the same time, to be sim¬ 
ple enough to find feasibility within current experimental 
technologies, for instance in circuit QED [53] or in sim¬ 
ulating all-optical setups [54|. Since non-Markovianity 
is linked to a dynamical recovery of the quantum coher¬ 
ence of a qubit mm, our work highlights that engineer¬ 
ing and exploiting suitably structured compound envi¬ 
ronments can supply useful developments for controlling 
and preserving quantum memory resources. It also mo¬ 
tivates further studies regarding the effects of multiple 
environments on the dynamics of correlations in many- 
qubit systems. 
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